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Tidal torque induced by orbital decay in compact object binaries 
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ABSTRACT 

As we observe in the moon-earth system, tidal interactions in binary systems can lead to angu- 
lar momentum exchange. The presence of viscosity is generally regarded as the condition for 
such transfer to happen. In this paper, we show how the orbital evolution can cause a persistent 
torque between the binary components, even for inviscid bodies. This preferentially occurs at 
the final stage of coalescence of compact binaries, when the orbit shrinks successively by 
gravitational waves and plunging on a timescale shorter than the viscous timescale. The total 
orbital energy transferred to the secondary by this torque is ~ 10^^ of its binding energy. We 
further show that this persistent torque induces a differentially rotating quadrupolar perturba- 
tion. Specializing to the case of a secondary neutron star, we find that this non-equilibrium 
state has an associated free energy of 10^^ — 10** erg, just prior to coalescence. This energy 
is likely stored in internal fluid motions, with a sizable amount of differential rotation. By 
tapping this free energy reservoir, a preexisting weak magnetic field could be amplified up to 
a strength of « 10^^ Gauss. Such a dynamically driven tidal torque can thus recycle an old 
neutron star into a magnetar, with possible observational consequences at merger 

Key words: binaries : close, gravitational waves, stars: neutron, stars: interiors, methods: 
analytical 



1 INTRODUCTION 

The effects of tidal interactions between celestial bodies have been 
known and studied for a very long time (Darwin 1879, Chan- 
drasekhar 1933, Kopal 1959, 1968). These effects range from the 
spin-orbit coupling - as observed in the earth-moon system- to the 
tidal disruption of astronomical objects, as they wonder too close 
to each other. 

In a binary, it is well known that a tidal torque can arise be- 
cause of viscous processes. In this picture, orbital angular momen- 
tum is transfeiTed to the stellar spin. If the process is efficient, the 
spin frequency may eventually equal the orbital frequency and the 
torque will vanish. This configuration is a caUed "tidally locked" or 
"synchronized" binary. This is typically observed in binaries con- 
taining normal stars (Zahn 1977). 

In two classical papers, Kochanek (1992) and Bildsten & Cut- 
ler (1992) showed that, in compact object binaries containing a neu- 
tron star (NS), there is not enough energy in the tides to lock a NS 
spin to its orbital motion, regardless of the magnitude of viscosity. 
Indeed, these authors further showed that all possible estimates for 
the NS viscosity would give only negligible spin energy. 

In this paper, we show that, when the inspiral is driven by 
gravitational waves, a torque of dynamical origin can ensue, even 
in absence of viscosity. The magnitude of this torque exceeds that 
due to viscosity in compact object binaries, at their final stage of 
coalescence. Because of the negligible effect of viscosity, the dif- 
ferential nature of the tidal torque induces a radially dependent ro- 



tating quadrupole perturbation. We quantify the total energy and 
angular momentum carried by this perturbation and its radial struc- 
ture. We find that it does not correspond to a minimum energy level. 
The free energy available is typically 10*^ — 10** erg, which may 
have interesting observable consequences. 

The paper is organized as follows. In Sections 2, we show that 
the dominant tidal torque prior to coalescence of two compact ob- 
jects is that produced by their orbital evolution. For regular stars, 
instead, the tidal torque is produced by viscous processes only. In 
§3 we derive the tidal radius and the radius where the dynamical 
instability sets in. In §4 we summarize the ordering of the charac- 
teristic radii, for different binary systems. In §5, we calculate the 
total amount of energy in the tides. In §6, we show the radial de- 
pendence of the rotating quadrupole excited in the secondary. In §7, 
we specialize to the case of a neutron star, described by a polytropic 
equation of state and we quantify the free energy stored in the ro- 
tating quadrupole. Finally, in §8, we discuss a possible outcome of 
this free energy and we draw our conclusions. 



2 THE INTEGRATED TIDAL TORQUE ONTO THE 
SECONDARY 

As it was shown in classical works dating back to Darwin (1879) 
and Chandrasekahr (1933), tidal interactions in binaries can ex- 
change angular momentum between the orbit of the system and 
each individual star. Let's denote Mi the mass of the primary star. 
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Figure 1. Sketch of the tidal torque which arises in a binary, when a tidal 
bulge is misaligned. The torque is visualized in the orbital plane, in the 
frame of the secondary (the larger, gray star). In this frame, the primary (the 
small black star) rotates about the secondary at the orbital frequency (we 
consider non-spinning star's). At the location of the secondary the dashed 
line indicates the contour of an unperturbed, spherical star. 



M* the mass of the secondary and i?* its radius. Tidal forces cause 
a departure from sphericity, which in the secondary's local frame 
corresponds to a non-resonant rotating quadrupolar tide (Thome 
1998). For brevity, we will frequently refer to this quandrupolar 
perturbation simply as the tidal bulge. 

In this paper, we work in the approximation of small tidal de- 
formations, which are valid for orbital separations a larger than the 
tidal radius 

At oc R*lq , 

where q is the mass ratio q = M,/Mi. The numerical coefficient 
is ~ 2 and it will be calculated as a function of the stellar internal 
structure in in Sec. 13. II 

If the tidal deformations are not perfectly aligned with the line 
joining the star centres (thereafter "line of centres" in short), each 
component will exert a non-zero torque on the companion and a 
flux of angular momentum is set up in the system. In the presence 
of a finite viscosity, there is a coupling between the orbital angu- 
lar momentum and the stellar spin. For compact object binaries, 
however, this viscous coupling is not efficient, and the angular mo- 
mentum exchange is between the orbit and the stellar bulge only 
(see section lXTt . 

We sketch here the classical derivation of the total tidal torque 
(A^t) applied onto the secondary, when its bulge is misaligned by 
an angle a (see Chandrasekhar 1933, Alexander 1973 and refer- 
ences therein). We assume small lag angles throughout the paper 
and, when needed, we will expand the trigonometric functions ac- 
cordingly, without explicitly mentioning it. We will show that the 
small angle approximation is justified for a ^ at and for typical 
stellar dynamical viscosity. 

We consider a two dimensional problem, in the orbital plane of 
the binary (see Fig.l). There, the primary at a distance a generates a 
tidal field -j/jt at the location of the secondary given by (Alexander 
1973 and references therein) 
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where the polar coordinates (r, 6) are centered in the centre of mass 
of the secondary, and 6 is measured from the semi-major axis of 
the bulge. In eq.[T] < r ^ i?*. 

The secondary reacts to the tidal pull readjusting to a new hy- 
drostatic configuration on a timescale i* — 2-k Jlo^. This mass re- 



distribution generates in turn a non-spherical perturbation -i/;, to the 
secondary's gravitational potential at each external point r > R*, 



V'* = «:2 ( ^ 
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(2) 



where K2 -the Love number- measures the global deformability 
of the secondary (Love 1909, Hinderer 2008). Because of ?/)», the 
secondary will apply a torque onto the primary given by Vetp* (r = 
a) X a. This is equal, but opposite in direction, to the tidal torque 
that the primary exerts onto the secondary. Therefore, 
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The value of a in the above expression must be regarded as an 
ejfective lag angle: a global property of a star with radius i?, and 
total deformability K2 . 

The torque vanishes when either the secondary or the sec- 
ondary's tidal bulge rotates with an angular frequency equal to the 
orbital frequency (i.e. a = 0). In the following sections, we dis- 
cuss the conditions for this synchronization not to happen and its 
physical consequences. 

2.1 The tidal interaction of a binary at & fixed separation 

2.1.1 the inviscid case 

Consider the simplest case first: a binary at a fixed separation a, 
and a non-spinning secondary. We work in a non rotating frame 
centered onto the secondary. 

In the absence of orbital motion, the tidal bulge is formed in 
a timescale i, and it lies perfectly aligned with the primary's po- 
sition. As soon as orbital motion is switched on, the primary be- 
gins to rotate about the secondary at the orbital frequency ujc, — 
y/GMtot/ofi, where Mtot = Mi + Af « . The bulge is initially at 
rest for a reaction time t*. In this time interval, the bulge lags be- 
hind the line of centres by a linearly increasing angle a ~ ujot. As 
a grows, a tidal torque (A^t) ensues that transfers angular momen- 
tum to the bulge. This process will last until the bulge acquires 
the angular momentum Jb = IhOJh = Ih^o, and it co-rotates 
with the primary. Here ojb and /b are the bulge rotation frequency 
and "ejfective" moment of inertia, respectively. Thus we can write 
NtU ~ Jb. 

Once this angular momentum has been transfered, the bulge 
rotates at ujh ~ ^o and it is perfectly aligned with the line of cen- 
tres. The torque vanishes again, having reached a stationary con- 
figuratiorQ From the above we can derive an expression for the 
effective moment of inertia of the bulge 



Jb = UNt/ujo = 67r^K2/(l + qf{ujo/uj*fM,Rl, 



(4) 



where we used the fact that a grows linearly in time during this 
transient phase. Eq.|4]shows that 7b is much smaller than the stellar 
moment of inertia /, « A/*i?^, given its strong dependence on the 
ratio {ijjo/^*) which is by definition <C 1. The angular momentum 



^ The angular momentum Jb is taken from the orbit through the action 
of the torque A'^t ■ This causes the orbit to loose a tiny amount of angular 
momentum, thus slightly shrinking and spinning up. However the orbit has 
an effective moment of inertia 7o 2> ^b (as we show below). It is thus 
clear that N'y / Io = dio <SC A'^t /lb = ^h- Therefore, the bulge frequency 
quickly reaches the orbital frequency and the torque drops accordingly to 
zero. This mechanism is self-regulating and does not produce any instability 
of the system. 
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that needs to be transferred to the tidal bulge for it to co-rotate is 
then Jb = /b^o oc K'hB?^u)o{'^o/ijJ*)^- 



2. 1.2 the ejfect of viscosity 

In the ideal case, where the secondary is perfectly inviscid, the 
circulation energy and the angular momentum of the bulge stay 
constant once the co-rotation is reached, and there will always be 
alignement of the bulge with the line of centres. In fact, the sec- 
ondary has a finite - though small - viscosity, which implies a con- 
tinuous transfer of angular momentum from the tidal bulge to the 
stellar spin. Thus, the tidal bulge is spun up by the primary torque 
but simultaneously slowed down by the viscous torque. The net re- 
sult is that the bulge cannot keep pace with the orbital motion of 
the primary and it lags behind the line of centers by an angle Qv. 
This, in turn, implies the existence of a persistent tidal torque. An- 
gular momentum is thus continuously transferred from the orbit to 
the bulge by the tidal torque and from the bulge to the secondary's 
spin by viscous stresses. At equilibrium the two transfer rates are 
equal: A'r ~ Jh/tv, where fv ~ Ri/v is the characteristic viscous 
timescale determined by the kinematic viscosity u. This condition 
determines the equilibrium value of the lag angle. 



tew ~ 1-2 yr 



-1 -8/3 



1 rad/s 



2-K UJo/{ui^tv). 



(5) 



where we used eq.[3]and eq.|4l and the fact that Ov is small to write 

CJb -^ ^o- 

The above derivation of Qv for weak viscous stresses is our 
alternative derivation to the classical one, which considers the re- 
sponse to the tidal field of each fluid element in the secondary star 
(cfr. Kopal 1968, Alexander 1973, Cutler & Bildsten 1992 and ref- 
erences therein). The fluid element will feel a time- varying external 
potential, the tidal potential of the primary, oscillating at twice the 
orbital frequency, 2ljo. This drives small radial displacements x{t) 
of the element from its equilibrium position. The fluid element is 
also subject to a restoring force, the self-gravity, which has a typ- 
ical oscillation frequency cj,. Finally, a weak friction, iv 3> t., 
is contrasting its oscillation. The equation of motion of each fluid 
element is thus that of a forced, (under)damped oscillator 

x{t) + 2t:;'^x{t)+ujlx{t) = 4e^^'^°'. (6) 

The solution is a displacement x{t) oc g»(2"oi-av) ^^jjich lags 
in phase with respect to the driving forcing oc e''^""* by an angle 
Ov OC LUo/[{ujf —uja)tv] ~ uio/{(^ttv). Thus, We recover our result 
(eq.H). 

If we start with the system at rest, the equilibrium value for the 
lag angle, a^ ~ cjotiag, is reached after a transient linear growth on 
a timescale tiag ~ t*(t«/tv) "C t.. This confirms our picture that, 
in the presence of viscosity, the system reaches first an equilibrium 
in which it does not perfectly co-rotate. In fact, it lags behind by a 
small but finite angle 

Qv = {oJo — '^s)ilag, (7) 

where ojg is the secondary's spin frequency. 

2.2 The effect of orbital shrinking 

In the previous section, we considered a fixed orbital separation. 
However, binaries may lose orbital energy and angular momen- 
tum by emission of gravitational waves. As a consequence, they 
will be brought towards the plunging phase which leads to coalesce 
through a succession of circular orbits on a timescale 



5/3 ' 



(8) 



where m. = A'L/Mq (Peters 1964). 

Indeed, for solar mass binaries which have a separation smaller than 
a ~ 3-Rq, tew is less than the age of the universe and merger 
can occur. However, for regular stars, this separation is smaller 
than their tidal radius, and tidal disruption inevitably occurs before 
GWs can become efficient in removing orbital angular momentum. 
Therefore, to investigate the effect of orbital decay by GW emis- 
sion, we should consider hard binaries, hosting compact objects. 
Their orbital shrinking changes substantially the classical picture 
of tidal interactions in binaries depicted earlier. 



2.2.1 the inviscid case 

We first ignore viscosity in the secondary star and consider, as a 
starting configuration, a binary in which the secondary's tidal bulge 
rotates exactly at the orbital frequency. As the orbit shrinks and, 
thus, the orbital frequency increases, the tidal bulge finds itself lag- 
ging behind and therefore subject to a tidal torque. This torque tries 
to spin it up to track the evolving orbital motion of the primary. 
As discussed above, the spin-up of the bulge occurs on a timescale 
~ t,. 

The lag angle induced by the orbital shrinking, cvgw, can be 
quantified in different ways. The easiest is to calculate the angle 
travelled by the primary on a time i*, as its orbital frequency grows 
at a rate cbc while the bulge lags behind at the original ljq: Qgw = 
(l/2)d;ot^ ^ LJ^/iuj^tGw). 

For an alternative, insightful derivation of ogw, we consider 
the dynamics of the system. The orbit loses angular momentum 
to GW emission, and we call the associated torque Nqw. As a 
consequence, the orbit spins up and a tidal torque A^t arises. This 
torque subtracts further orbital angular momentum, transferring it 
to the secondary's bulge. Equilibrium is reached when the bulge 
frequency becomes able to change at the same rate as the orbital 
frequency, ojo = ujb- This requires the existence of a non-zero tidal 
torque - and thus of a non-zero lag angle qgw-, which continu- 
ously spins up the bulge at the appropriate rate, 

Nt _ Now + Nt 



/o 



(9) 



where Jo is the binary moment of inertia. Since the orbital decay is 
largely driven by GW emission, we neglect the second term on the 
right-hand side of equation |9] and solve for the equilibrium angle 



QGW — 27r 



^*iGW 



(10) 



Finally, in analogy with section 12. 1. 21 we can derive qgw 
from the oscillatory motions of fluid elements in the secondary 
star, which are forced at an increasing orbital frequency. Differ- 
ently from eq. |6] we neglect here the damping term but allow for a 
time-dependent forcing term, 

x{t) + ujlx{t) = Ce ^ '"^'''•'^''' . (11) 

The solution for the radial displacement is x(t) oc e''-'^ 2u;a(t)dt^a] 
Therefore, even in the absence of viscosity the evolving orbital fre- 
quency will cause the displacement of each fluid element to lag 
in phase with respect to the perturbing potential by the angle a. 
To first order in (ajo(i)/<^»)^, we get a = qgw oc luo/[{uj^ — 

i^o)tew] ~ Ulo/{ujftGw). 

From these arguments, we conclude that, even in a perfectly 
inviscid secondary, the tidal bulge will develop an angular lag and 
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will thus be continuously subject to the minimal tidal torque asso- 
ciated to it (Lai, Rasio, Shapiro 1994b, Lai & Shapiro 1995). 



2.2.2 In the presence of viscosity 

In the presence of viscosity, there is a net transfer of angular mo- 
mentum from the tidal wave to the secondary's spin. In section 
12. 1.2l we showed that, for a given ujo, the value of Ov is determined 
by the balance between the viscous torque, draining angular mo- 
mentum from the rotating bulge, and the tidal torque, that acts to 
replenish it. If there is orbital shrinking, but it proceeds at a slow 
pace, i.e. taw S> iv, the effect of ojo can be neglected and the 
picture of section l2.1.2l is still the valid description of the equilib- 
rium configuration. If, on the other hand, icw ^ iv, the lag angle 
is instead determined by the orbital evolution as described in the 
previous section. Note that tow <^ tv implies Qv <C acw and 
viceversa (compare eq. |5]and eqllOt. This means that during the 
binary evolution the effective tidal torque is the maximum between 
the "viscosity-induced" and the "GW-induced" torque. 

Since tew oc a*, we expect the GW-induced torque to be 
important at small binary separation. Indeed, acw > "v when 



u)o > 1.1 ni-^ 



5/8 ( tv_ 

yr 
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(12) 



where mi = Mi/Mq. In ordinary stars, the effective viscosity 
operates on a timescale of ~ yr (Zahn 1977, Hut 1981, Goldreich 
& Nicholson 1989). Therefore, an equal mass binary with a so- 
lar type companion will never enter the high frequency regime of 
eq. [T2] since it requires a < 6 x 10 cm <^ Rq. The estimated 
viscous timescale in a typical NS ranges from tv ~ 0.1 s for the 
crust anomalous viscosity (Kochanek 1992), to tv ~ 1 yr for shear 
viscosity (cfr. Andersson 2007). For WDs the shear viscosity acts 
on a very long timescale fv ~ 10^^ -^ 10^** s (cfr. Neubauer 1986). 
The internal magnetic field has been proposed as a possible source 
of anomalously high shear viscosity in WDs, giving a somewhat 
shorter t^ « 10® s. These ranges of timescales together with the 
large compactness of these objects allow for the condition eq. [12] 
to be satisfied for binary separations much larger than their tidal 
radius. This can be easily seen by recasting eq.[T2]in terms of a, 

"■ ,„4 -2/3 1/3 / 12Km W f, y/" 

< 10 m^ mj \ — ;:; — I I TTT I ' (13) 
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3 TIDALLY-INDUCED DYNAMICAL INSTABILITY 

When the orbital separation in a binary becomes comparable to the 
size of its components, tidal effects cause a non-negligible steep- 
ening of the interaction potential with respect to the term oc 1/a. 
This increases the radial component of the gravitational force, and 
alters the keplerian relation between orbital separation and period. 
Eventually, a point can be reached where no closed orbits are pos- 
sible and the system becomes dynamically unstable. The binary is 
brought to coalescence in a dynamical (free-fall) time, much faster 
than dictated by GW emission (Lai & Shapiro 1995 and references 
therein). The same arguments that lead to the definition of acw 
thus lead to the definition of a new lag angle, Qdyn. This is obtained 
from the former, by simply substituting icw with the dynamical 
time in the relevant regime. 

In order to determine when the above instability occurs, it is 
first necessary to calculate the star deformation as a function of the 
orbital separation. In practise, we calculate the evolution of 77 = 



h/ Rt, the fractional hight of the tidal bulge, which allows us to 
derive the perturbed stellar radius, 

R{a) = R,[l + Ti{a)]. 

In fact, this is the semi-major axis of the star: the elongation along 
the line of centres (9 = 0). To derive r), we will impose (in Sec. 3.1) 
that at each radius the secondary's total energy is equal to that for 
a point mass (the tidal term causes just a redistribution of energy). 
This calculation will also allow us to find an analytical formula for 
the tidal radius ot and determine the radius, Qc, at which the binary 
comes into physical contact. Then (Sections 3.2 and 3.3), we will 
write the total (newtonian) effective potential for the binary system 
and find at which separation a minimum cannot be found anymore. 
This separation ad marks the onset of the dynamical instability. In 
all cases, our conclusions compare well with published results (cfr. 
Clark & Eardley (1977), Lai & Shapiro (1995), Bini, Damour & 
Faye (2012) and references therein). 



3.1 Tidal stellar deformation and tidal radius 

The secondary's tidal bulge is a result of the work done by the 
first term in the expansion of the primary's newtonian potential 
ipT (eq[T). This causes a decrease in the binding energy of the 
originally unperturbed secondary. In turn, the change in the sec- 
ondary's external potential, 7/;* (eq.O, is a result of the tidal bulge; 
it increases the binding energy of the orbit, further drawing from 
the stellar binding energy. In the system's energy budget, the addi- 
tional energy terms t/jt and ijj-t are thus balanced by the decrease in 
the stellar binding energies. This statement of energy conservation 
reads, 



KU^T + Mi?/), 



PGM'i I3GMI 



R* 



R 



(14) 



where tpT and 1/;* are evaluated along the line of centres (9 = 0), 
for a separation r = a and a stellar radius R (instead of i?,). The 
structural constant in the self gravity term (right-hand site) is /3 = 
3/(5 — n) for a polytropqj of index ?i (cfr. Lai & Shapiro 1995). 

A manipulation of eq. [14] leads to the following algebraic 
equation for rj. 



K2(l + r?)® + g 



(l+vf-q^P 



77 = 0, 



(15) 



For a given polytropic index n and the mass ratio q, we can solve 
numerically the above equation to obtain i]{a). In Fig.|2]we show 
our results for a NS-NS system, adopting different values of q and 
n. In the left panel, we show a binary system where both NSs are 
described by an n = 1 polytropic equation of state (eos) , a good 
approximation for fiducial NSs with M=1.4 Mq and R=12 Km (cfr. 
Sec.|7). In the right panel, we consider a primary NS which is sig- 
nificantly more massive, ~ 2 Mq . This case can be represented by 
a stiffer eos with n = 1/2 polytropqj 

From the above results, we obtain an unambiguous definition 
of the tidal radius as the orbital separation below which the residual 
binding energy of the secondary is unable to balance the tidal field. 
This always occurs when 77 reaches a maximum value of r]T ~ 
0.4, independently of q and barely sensitive to n (see Fig. 2). This 
deformation corresponds to a ratio of semi-minor, R^,{l — 7?/2), to 



^ In SecQwe will discuss polytropic models in greater detail. 

•^ In a BH-NS system only the secondary, i.e. the NS, would be deformed. 

This is just a sub-case of the ones plotted in Fig.|2] 
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Figure 2. Fractional tidal elongation, r) fs h/Rt , of the radii of both components in a NS-NS system, as a function of the orbital separation. Left Panel: two 
NSs described by the same e.o.s., for three different values of their mass ratio. The black lines are for the secondary, the blue ones for the primary star. Note 
that for g = 1 the two curves are exactly superimposed, since the NS are identical. Right Panel: the secondary NS is described by an n = 1 polytrope and the 
primary one by a stiffer, n = 1/2 eos. Note that, for g = 1, the stiffer NS has a larger deformation, at a given a, and is disrupted at a slighlty larger orbital 
separation. This is because it is less centrally condensed, and its Love number is thus larger. On the other hand, for q ^ 0.95 the less massive star is more 
deformed and tidaly disrupted, dispate the steeper density profile. 



semi-major, _R, (I + 77), axis of 0.55. Indeed, no real solution toeq. 
[T5]exists below the minimal orbital separation ax corresponding to 



(2.14 



2-30);^. 



(16) 



where the numerical coefficient depends only on n, with the small- 
est coefficient obtained forn = 1 and the largest one for a n = 1/2. 
We can thus confirm the classical scaling, ax oc q~^'^. Note that 
this expression for the tidal radius holds in general, whether the 
primary is a NS or a BH. 

Given the simplicity of our treatment, eq. [16] compares re- 
markably well to recently published results of relativistic numer- 
ical simulations (Taniguchi et al. 2008, Pannarale et al. 2011). For 
example, Pannarale et al. (2011) find that the tidal radius in their 
simulations corresponds to a stellar axises ratio of 0.44, which is 
not so different from our newtonian result (0.55). Indeed, our def- 
initions (eq. [14] and list provide a simple physical argument for 
interpreting published numerical results. 

Finally, we conclude this section by comparing our results for 
»7 with the classical approximation h/R, = tf}T/{l3GM,/R,), 
where <^* is not included (solid red line in Fig. 2). This approxi- 
mation clearly fails to describe the final stages of the inspiral and 
to determine the value of the tidal radius. 



3.2 NS-NS binaries 

We can now describe the final stages of coalescence of two NSs, 
tracking the growth in radius of both components. For convenience, 
we will label Ri and R2 the perturbed stellar radius of the primary 
and secondary star, respectively. 

Beside the tidal radius (ea.ll6t, there is another characteristic 
radius at which the inspiral may end: the orbital separation Oc at 
which the NSs first come into contact. The question is which one 
is the larger. In Fig. [3]we show ht and ttc as a function of q for a 
specific case of a binary with identical NSs. For mass ratio q > 0.7, 
ic <; clt- In fact, this is a quite general result. At the tidal radius 
77 ~ 0.4 for the secondary star, implying that its perturbed radius 



is i?2 ~ lAR,. The primary's radius will also be deformed, by an 
amount which is ~ q'^ times smaller than the secondary's. The sum 
of the two radii will thus be i?2(aT) + i?i(aT) ~ (2.4 + 0.4(7^)i?*, 
which for g > 0.7 is larger than the tidal radius for any polytropic 
index. For smaller values of q, the tidal radius would instead be 
reached before contact. However, mass ratios smaller than 0.7 are 
highly unlikey for NS binaries, and we do not consider them here. 
We thus conclude that physical collision will generally occur prior 
to tidal disruption in NS-NS systems. 

The last characteristic separation to compute is ad, where the 
system becomes dynamically unstable. The total effective potential 
for two NSs is 



L'^jl + q) _ ALMi 
2M.a2 a 



'-^^iRl + Rl) 



.(1) 



M^Rl 



(17) 

a^ a" 

Here the first term is the centrifugal potential, the only repulsive 
term in the sum, where L'^ — L2 + L\ is the total angular momen- 
tum. The second represents the interaction potential between point 
masses. The third and forth are the two leading-order terms in their 
mutual tidal interaction (i/jt and ip*), evaluated along the line of 
centres (6 = 0). In the above, the perturbed stellar radii R\,2 are 
a function of a (see previous section), but, we omitted to write the 
a-dependence for simplicity. 

The above potential will have a minimum - coiTesponding to 
a stable, circular orbit - if and only if there exists a point where its 
first radial derivative is zero and its second derivative is positive. In 
the opposite case, there will be no minimum and the system will 
be dynamically unstable (cfr. Clark & Eardley 1977). Formally, 
derivatives of R with respect to a should also be taken to derive 
the instability limit in this way. Such a rigorous approach would 
greatly increase the difficulty of the calculation, adding little physi- 
cal value. In fact, it is possible to show that for 77 ^ a few tenth, the 
terms with dR/da can indeed be neglected. |j This approximation 

■* Let's take the derivative with respect to a of any of the terms in ~ R? /a^ 
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will lead us to find that the dynamical instability sets in for 77 ^ 0.1, 
justifying a posteriori our procedure. 

With our approximation well defined, we now impose the first 
derivative of eq.[T7]to be zero, 






M,Mi 



+ 



a 



+ 



GK'^'MiRl 



3MiM. 



[Rl + Rl), 



a' a* 

while requiring the second derivative to be positive, 

3L2(l + g) 



(18) 
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42K,i^'>MlR^ 



42Ki^^M^Rl 
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UMiA'L 



(Rl + Rl). (19) 



Combining the two conditions above and after some manipulation, 
we obtain an equation for the adimensional radial distance, j/d = 
a-d/R2{a.d), where the instability sets in: 



yl-3yl 



' + || 



24 



..(2) 



\^ 2 4^ fRl 



.(20) 



To test the validity of our approximation, we considered the 
specific case studied by Lai & Shapiro (1995) of a NS-NS binary, 
where the stars have the same radius Ro = 12 km, polytropic index 
n = 1/2 and have a mass ratio q = 1/2. These authors solve nu- 
merically the full set of (newtonian) hydrodynamical equations and 
find that the instability ensues when a^/R* — 3.1 (cfr. top panel of 
their Fig. 2). Our analytical method (eg. 120) gives ad/R* — 3.05. 
Considering the significantly different procedures, the agreement is 
quite remarkable. 

We can now proceed to solve equation |20] in other cases of 
interest. In Fig[3] we plot the solution as a function of the mass 
ratio, when the secondary is described by an n = 1 polytropic eos. 
It shows that the dynamical instability always occurs before the 
two NSs come into contact (cfr. Lai, Rasio & Shapiro 1994a; Lai 
& Shapiro 1995). 



3.3 BH-NS binary 

Before repeating the same calculation for a BH-NS binary, we 
briefly address the notion of a "contact" radius in such systems. We 
consider the event horizon (_Rh = 2GMi/c^) as the relevant radius 
for a (non-spinning) BH, and assume that the secondary will disap- 
pear in the event horizon once its elongated radius first touches it. 
Thus, the contact radius Oc in a BH-NS binary is the orbital sepa- 
ration at which the NS is swallowed in the event horizon. 

We now ask what is the condition for the tidal radius to be 
reached before contact, so that tidal disruption will occur in an in- 
teresting region. The sum of the two radii when the tidal radius ax 
is reached can be written as 



R2{aT) + Rh ~2ii. 



0.7 +- 



(21) 



in eq. [IT] We then impose 2{R/a^)\dR/da\ -S (SR^/a"') or equiv- 
alently |d?7/da| <g (3/2) (1 + r])/a. For a power-law dependence of 
rj oc {Ri,/a) , this conditions becomes r] <^ 3/(25 — 3). For any rea- 
sonable 5 (5 < 10) this implies r; <C a few tenth. A consistent result is 
obtained considering the terms in R^ /a® . 
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Figure 3. Tidal radius (ax ), contact radius (Oc ) and radius of dynamical 
instability (a^) as a function of the mass ratio q, for a NS-NS binary where 
both components are described by an n = 1 polytropic eos. Dynamical 
instability always sets in before contact. The specific eos was chosen for 
the sake of illustration, but our conclusions have a general validity for any 



where we defined the NS compactness C — GMt/{c^Rt) and 
used the fact that tjt ~ 0.4. Using eq. J16t . and after some manip- 
ulation, the condition ht > flc leads to the inequality. 



q- (1.53^1.64)5 



2/3 



+ (&^°- 



(22) 



The above condition can be fulfilled for g-values that depend 
on the NS compactness. For our fiducial NS, the relevant coeffi- 
cient is the smaller one. This gives q > 0.11, or Mi < 12.7 Mq. 
On the other hand, the fiducial NS will disappear in the BH hori- 
zon before being disrupted if the BH is more massive than ~ 12.7 
Mfj. Changing the NS structure gives slightly different numerical 
values, without affecting our conclusion that for stellar mass black 
holes, the tidal disruption generally occurs outside the event hori- 
zon. 

A well-known property of the Schwarzschild potential is its 
steepening at close orbital separations. This genuinely relativis- 
tic effect determines the so-called innermost stable circular orbit 
(ISCO). Although different in nature, this is totally equivalent to 
the dynamical instability discussed before, for our purposes. In- 
deed, once this separation is reached, orbits cannot be closed any- 
more and the secondary will plunge into the event horizon on a 
dynamical timescale. Therefore, we also want to establish the lo- 
cation of the ISCO with respect to the tidal radius, in order to de- 
termine its role in the orbital evolution of a BH-NS binary. For a 
Schwarzschild BH, aisco ~ SRh- Using eq.[T6]for the tidal radius 
it is straightforward to see that the ISCO will be larger than the 
tidal radius if g < (0.3 + 0.33) (C/0.17)^/^ for 1/2 s^ n sC 1. 

Finally, we discuss the condition for dynamically instability 
in BH-NS systems. The reasoning is identical to that in the previ- 
ous section but the formulae must be slightly modified to account 
for the fact that the BH does not suffer tidal deformations (its Love 
number is zero). As a consequence, we expect the dynamical in- 
stability to set in at smaller orbital separations, for a given set of 
system parameters. It is a simple matter to re-write eq. |20]for a 
BH-NS system and obtain. 



5 o 3 o/l '^2 

Vd - 3yd = 24 — , 

q 



(23) 



The solution of this equation as a function of q is plotted in 
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Figure 4. Tidal radius (ax), contact radius (ac), radius of dynamical instability (a^) and ISCO as a function of the mass ratio q, for BH-NS binaries. Left 
Panel: the secondary NS has M, = 1.4 Mq, R, = 12 km and n = 1. Right Panel: the secondary NS is a massive, 2 Mq NS with R, = 12 Km and an 
n = 1/2 eos. 



Fig.m along with the tidal radius, the contact radius and the ISCO, 
for two different BH companions. The main result is that for a BH 
with masses greater than a few Mq , the relativistic dynamical in- 
stability sets in before the tidal one (i.e. fliaco > id)- In this regime, 
the secondary is generally tidally disrupted -as opposed to swal- 
lowed - for stellar mass black holes. 



4 SUMMARY OF THE EVOLUTIONARY SCENARIOS 

In the previous sections we analyzed several physical mechanisms 
that produce important non-trivial effects in the late evolution of co- 
alescing binaries. These become important at different orbital sep- 
arations. In the following, we summarize our main results for the 
ordering of these phases, in NS-NS and NS-BH systems. 



4.1 NS-NS binaries 

The situation is relatively simple in these systems. The orbit be- 
comes dynamically unstable at the orbital separation ad, which is 
larger than any other relevant size. Subsequently, the two NSs will 
inevitably touch before reaching the tidal radius (eg. 116). This gives 
the following ordering: ax < ac < ad. 

The orbital evolution of a NS-NS system will thus be driven by 
the emission of GWs down to Od. Once in the dynamically unstable 
region, the binary orbital shrinkage is significantly accelerated and 
occurs on the dynamical timescale (id) at as_. The NSs eventually 
touch at ac. The inspiral ends here and the complex merger phase 
sets in. Correspondingly, we expect the tidal torque to be charac- 
terized by a lag angle acw during the first part, and by a larger lag 
Qd ~ Q!Gw(iGw/id) once the dynamical instability sets in. This 
will be discussed in greater detail in the next sections. 



4.2 BH-NS binaries 

The mass ratio in these binaries turns out to have a large effect 
in determining the relative locations of the characteristic radii. The 
dynamical instability occurs in these systems only for relatively low 
mass BHs. However, we noticed that the ISCO can play a similar 
role at significantly larger orbital separations, when the BH mass 
grows. The contact radius is generally smaller than the tidal radius. 



and only for very small mass ratios (< 0.1) it plays a role. We 
neglect such small q values here. As a consequence, we expect three 
types of orbital evolution for BH-NS binaries, depending on the 
mass ratio. For a fiducial NS mass of 1.4 M0 and radius f 12 Km, 
this will be a function of the BH mass only. 

• if Ml < 4 M0 we have the ordering aisco < ax < ad. In 
this case, GWs drive orbital shrinking until the dynamical instabil- 
ity is reached. A short plunging phase follows, until the NS reaches 
the slightly smaller tidal radius at which it gets tidally disrupted. 
The orbital evolution stops at ax; 

• if Ml > 4 Mq the ISCO will coiTespond to the largest orbital 
separation. In this case, GWs drive the orbital evolution until the 
NS reaches the ISCO. A plunging phase follows also in this case, 
whose radial extent grows quickly with increasing Mi. The orbital 
evolution will again be stopped at ax ; 

• if Ml >14 Mq we get ad < ax < ac < aisco- Also in this 
case, GW-driven orbital shrinking will be overtaken at aiaco by a 
plunging phase. The NS is swallowed in the horizon before being 
tidally disrupted. 



5 THE TOTAL TIDAL ENERGY 

The bulge has angular momentum but also a rotational energy given 
by E^a = (l/2)/baj^. If not dissipated, this energy remains stored 
in the wave and increases as the orbit shrinks. We now evaluate 
the total amount of energy Eh.T that the tidal torque transfers to 
the secondary's bulge, as the system evolves towards coalescence. 
Since the lag angles considered are small, we approximate ujh — 
ujo in order to simplify the calculations. This implies that the gain of 
energy is diSb = Hh^oduio, regardless of the speed of the inspiral. 
Thus the total energy released up to an orbital frequency ojq is 



-Bb,T(^o) = 



3TV^K2M,R'i 

{l + qy^ji 



uj^dujo 



(24) 



where we assumed an infinite initial orbital separation. This is justi- 
fied by the fact that essentially all of the energy is transferred close 
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to the final radius, and thus the starting point of the integration is 
irrelevant. In the most interesting cases, the inspiral ends at the tidal 
radius, where the orbital frequency is aj(aT) = ijJT, 



UJT = ci a;*(l + qy'^ ^ (2.7 -^ 3) X 10' 



M. 



1.4 M 



1/2 



(25) 



12 Km 
R, 



3/2 



(1+g)^/^^ 



In the above, the coefficient ci = (0.29 -f- 0.32), where the smaller 
value is for n — 1/2 and the larger is for n = 1. Thus, from eg. 1241 



£b.' 



; 3 X 10" 



E, 



K2 



4.3 X 1053 0.26 






erg, 



(26) 



where _B* = GM^ / R, is the secondary binding energy. The Love 
number hi2 was normalized to a typical value for an n = 1 newto- 
nian polytropqj and the orbital frequency at the tidal radius is given 
by eq.|25]with c\ — 0.32 for an ?i = 1 polytropic NS. 

The total energy transferred to the bulge in an equal mass bi- 
nary is thus ^ 10^3 — 10^^ of the secondary binding energy. The 
numerical coefficient in the last line of eq.|26]is for our fiducial NS 
parameters. A solar mass WD, with _R, = 5 x 10"^ km, has instead 
E,Ki5x 10™ erg and £b,T ^ 1.3 x lO'*** erg (for WDs K2 - 0.1, 
cfr. Verbunt & Hut 1983). 

The orbital energy i5b,T(i^T) which is transferred to the sec- 
ondary can be of the same order of that required for co-rotation at 
at: ~ (l/2)(0.35)A/.i?^a;^ ^ 0.02E,. However, Et^^T is mostly 
stored in the quadrupolar oscillations within the secondary and only 
a minor fraction (iow/iv) is actually viscously converted into spin 
energy. We, therefore, confirm that binaries with compact objects 
cannot reach co-rotation before they merge (Kochanek 1992; Bild- 
sten & Cutler 1992). 

In the following, we will show that this internal motion, i.e. 
the rotation of the tidal bulge, is a function of radius within the 
secondary. In this respect, it inherits the differential nature of the 
tidal force. 



6 THE DIFFERENTIAL STRUCTURE OF THE TIDAL 
BULGE 

Our description of the tidal bulge is that of a wave which rotates 
in the secondary's frame at the ratqjcjb ^ i^o, exciting quadrupo- 
lar fluid oscillations. As such, it carries a total angular momen- 
tum, Jb (§ 12. It and rotational energy Eh (§ HJ. We can thus as- 
sociate to the wave an "effective moment of inertia" /b (eq. ^, 
such that Jb = /bij-ib and Eh = {l/2)Ih'-^h- If circumstances 
arise that lead the wave to rotate slower than the tidal potential, 
the system seeks immediately a new equilibrium configuration, ac- 
celerating the wave, via a tidal torque A'^t. We showed that this 
happens when GWs drive orbital acceleration iJJo, or when the orbit 
becomes dynamically unstable, at small orbital separations (either 
at Qd or ftisco). In those cases, we argued that mechanical equilib- 
rium requires the tidal torque to accelerate the bulge at the same 
rate, Uh ~ iJJo- We now proceed to calculate the radial structure of 
this rotating quadrupolar tide (or wave). To this aim, we divide it in 



5 Relativistic models give 0.05 J- K2 < 0.1 (Hinderer 2008; Postniakov 
2009). We adopt the newtonian value to compare with our results of ijT] 
^ Strictly speaking, it is ~ oJo — i.Js, but we are neglecting the secondary's 
spin. 



concentric rings and assign to each ring the mass contained in the 
corresponding spherical shell of the NS. We then require each ring 
to be in mechanical equilibrium. Given the marginal role played 
by viscosity in the systems of interest, we can assume that there 
is no transfer of energy and momentum both among the rings and 
between the tidal bulge as a whole and the stellar spin. 

The total tidal torque A'^t depends on 3 global quantities of the 
secondary (see eq|3}: its radius, Love number K2 and lag angle a. 
Hence, the same expression holds for any radius r < R,, with the 
appropriate global values of K2 and a. The difference between the 
torques applied to two concentric portions of the star, with radii r 
and r + dr, will thus constitute the torque applied to the ring with 
inner radius r. To first order in dr, we get 



dNrir) = 3^^^d{K2r''a) 



(27) 



where we leave indicated a generic lag angle a. 

The equilibrium condition demands that this torque acceler- 
ates the ring so to track the orbital acceleration: ijh{r) = uJo = 
const. The ring angular acceleration is ojb ~ dNT/dIt, where the 
right-hand term can be written as the ratio between 
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(28) 



(29) 



where the last equation is obtained by taking the derivative with 
respect to r of eq.|4] Therefore, 



^b 



dNr 



2n^Gf{r) = 6j^, 



(30) 



where /(r) is the ratio of the two radial derivatives in eg. 128 1 and 
eq.|29l Solving eq.|30]for the lag angle a, we obtain 

a{r) = 2-K^Uo/u}lir). (31) 

Generally, eq. 13 II implies a tidally distorted structure, where 
the lag angle varies within the star. This is because the restoring 
force, ujt{r) — GAI{r)/r^, is usually a function of r, where 
M{r) is the mass within r and M{R,) — M,. This conclusion 
is sketched in Fig. |5] A notable case is that of a constant density 
profile, which will ensure no distortion and a lag angle given by the 
global value a = 2-K^Lbo/uJt (eq.llOll. 

The exact dependence of a on r is determined by the func- 
tional form of AI{r). In order to proceed, we thus need to consider 
a specific stellar structure. 



7 TIDAL EFFECTS IN NEUTRON STARS WITH 
POLYTROPIC EOS 

In this section, we calculate the radial profile of a rotating 
quadrupolar tide, which is induced in a NS by the tidal field of 
the primary. 

For this purpose, we approximate NSs as polytropic stellar 
models, where the pressure P is given by P = Kpp^'^^'". Such a 
model, though idealized, describes reasonably well the NS global 
properties obtained with more sophisticated equations of state (Lat- 
timer & Prakash 2001; Lattimer & Prakash 2006). The equation of 
state, together with the stellar radius and mass, completely deter- 
mine the stellar profile. In particular, a n = 1 index is well suited 
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Figure 5. The structure of the rotating tidal bulge, resulting from the radial 
variation of the lag angle. The crest of the wave at different radii is in phase 
with different positions of the primary along the orbit. 



to describe a typical NSs withi a radius R, = 12 Km and mass 
M* = IAMq. However, recent indications for a 2 Mq NS (Lat- 
timer & Prakasli 2010 and references therein) suggest that a NS eos 
might range in stiffness. More massive objects would thus be bet- 
ter approximated by stiffer polytropes, with n = 1/2 or n — 2/3. 
According to the detailed analysis by Lattimer & Prakash (2006), 
most eos's compatible with current observational constraints give 
radii in a fairly narrow range for NSs in the (1 — 2) Mq range. 
Therefore, we adopt the reference value of i?, = 12 Km for the ra- 
dius for any NS. On the other hand, we associate n = 1 polytropes 
to 1.4 Mq NSs and either n = 1/2 orn = 2/3 to 2 Mq NSs. 



7.1 The differential rotation of the tidal wave 

For a polytropic star the equation of hydrostatic equilibrium is 
reduced to the Lane-Emden equation between the dimensionless 
radial coordinate ^ and the dimensionless density 6. ^ is related 
to the radius by r = A£^, where the constant A^ = (n + 



l)KJ{A7vG)p, 



(l-n)/n 



The dimensionless density distribution is 



6" = p/pc where pc is the central density. The outer boundary 
(Cout) is found as the first zero of the solution 6(^) of the Lane- 
Emden equation. Knowing ^out and the star's physical radius, one 
can derive A — R, /^out • 

Integrating 9 over ^ one obtains an expression for the mass 
within radius ^ as 



M(< = 47rylVc / C"e"(C')dC' = -47rAVcr 
Jo 



de 



(32) 



The characteristic frequency of oscillation of matter within radius 
^ can now be written as 






IdB 



de/dc 



C [dO/d^L 



(33) 



From this equation, we can obtain the corresponding value of the 
lag angle ogw when GWs drive the orbital evolution, 



aGw(C) = 27r 



2^o 



2-R 



2i^o F„ (^out ) 



(34) 



where -Fn(C) = i^/OdS/dS,. The run of (xJ,(r) from the centre to 
the outer boundary is shown in Fig.|6]left panel. The coiTespond- 
ing run of qgw(t') at a = ad is shown in Fig. [6]right panel. The 
external layers of the bulge react more slowly to the change in po- 
sition of the primary and lag behind the NS core, which instead, is 
promptly realigned. 

With an identical reasoning, we can define a generic dynam- 
ical lag Od, that holds when the orbit evolves on a dynamical 
timescale id. 



a4^)^2Tv' 



LoUOU 



(35) 



In Fig.[6]right panel, we plot ad at a = at (blue line). The absolute 
value and especially its derivative are larger than those of acw at 
Qd (black line): during the dynamical plunging the star becomes 
significantly more differentially perturbed. 

From the above expressions, we can define the rotation fre- 
quency of each portion of the tidal bulge, in analogy with the clas- 
sical description of the equilibrium tide in viscous bodies (eq|7]l. 
The lag time with which the bulge reacts is now i* (r), while ujh{r) 
plays the role of the spin frequency ajg. At a given radius r, the 
corresponding relation between acw (or Od) and a;b is thus 



[ujo - ujb]t,{r) = aGw{r) if a > ad 

[ojo — tJb] t, (r) — Qd('") if aT,ac^a<ad- 



(36) 



When the orbit becomes dynamically unstable, either at ad or at 
aisco, the system evolution will depend on the details of the plung- 
ing phase, whose treatment is beyond our scope. Here, we sim- 
ply assume that, at the onset of dynamical instability, the relevant 
time will switch from the GW-driven inspiral time (eq. |8) to the 
much shorter dynamical time calculated at ad : td ~ 2-k/u}^, where 
ujd = ct)o(ad). Once this occurs, our most conservative assumption 
will be to freeze the further growth of the tidal lag, fixing a^ to the 
value it has when the dynamical instability first sets in. With this 
assumption, the radial behaviour of ujh along the orbital evolution 
is completely defined. Inserting eqs.|33]|34]and|35]into eq.|36] gives 



^b(0 = 
^b(0 = 
wherqj 



^o [1 - ecw G'„(^)] if a > ad 
^o [1 - Ed G„(^)] if aT < a < ad 



(37) 
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2u!tU!o 



^n VSoutJ, 



and 



Gn(0 



F'J\i) 



-,1/2 



d6l/dC 



Eq. [37] shows that the tidal wave rotates in the secondary's frame 
at a non-uniform angular speed, which decreases from the center 
outwards (Fig. IT}. 



Note that ed/^GW = tGw/Po 2> 1, Po being the orbital period. 
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Figure 6. Radial structure for n = 1 polytropic equation of state. Left panel: the characteiistic frequency Ut ; Right panel: the lag angles oqw (black line) 
and Q(j (blue line) at respectively the dynamical and tidal radius. All quantities are plotted as a function of the normalized radius § = 7rr/i?, . The binary 
mass fraction is tj = 1. 



7.2 The total free energy 

Our results outline a new picture, which may have important con- 
sequences. As a result of the coupling between tidal interaction and 
orbital evolution, a radially-dependent tidal perturbation will be ap- 
plied to the NS. Fluid elements at different radii within the NS react 
with different time lags at the passage of the perturbing primary. 
The element at the surface, being the slowest to react, will still be 
pointing towards an earlier position of the primary, while a central 
element will be instantaneously aligned with its current position. 
In other words, oscillations of fluid elements at different radii are 
out of phase and occur along different directions. This condition is 
not consistent with hydrostatic equilibrium and it is not a minimum 
energy configuration. 

We thus anticipate that internal fluid motions will be excited, 
with both radial and azimuthal components of velocity (cfr. Goldre- 
ich & Nicholson 1989): these internal motions will carry that part 
of energy of the tidal perturbation which is in excess with respect to 
the equilibrium state. A complete characterization of such motions 
is beyond the scope of this paper. However, we can calculate the 
amount of free energy associated with them. 

To this aim, let us define some proper weighted average of the 
wave angular velocity, ujh- The frame rotating at this frequency will 
be in phase with the tidal wave at a particular radius. As a conse- 
quence, an observer in that frame would see elsewhere oscillations 
with a 0-component of the displacement. The amplitude of such 
^-component would grow with distance from the reference radius. 
The energy associated with the motion observed in this frame repre- 
sents a genuine free energy reservoir carried by the wave. If, indeed, 
such motions were to be damped by some, yet unspecified, internal 
torque, all rings would oscillate in phase and the wave would rotate 
at the (uniform) angular velocity Uh- The difference - at constant 
total angular momentum - between the energy in the excited state 
and the minimum energy corresponding to the uniformly rotating 
wave is the free energy reservoir that we want to quantify. 

The total angular momentum and energy caixied by the wave 
(the tidal bulge) can be expressed as 



dh 



iOt{S,)d^ = Ib.oi^oih — e/2), 



Eu 



1 f dib 2/fN ,<^ h,o^o , r 



where we used eq.|29]to obtain 



dh _ ,iidlb _ d 

1^' iR-^'^-'di 



FniO 



i-2eJ2 + £ /a), (38) 



(39) 



where /b,o = GTv^A'^u}i/\G'^(l + Q]'^4-KPr]. The integrals in eg. [38] 
are 



h 



d_ 
di 
d_ 

d! 

d_ 
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H2C 



K2e' 



.FniO 



FniO 



d^ = h/hfi, 



(40) 



The energy of a uniformely rotating wave is instead Emin ~ 
(l/2)Jba)b — itot/(2/b). Therefore, the free energy at a given 
orbital frequency is 



r 2 

dljjf = iitot — iimin — e 



r ^ 

13, :- 



(41) 



where e can be either ecw or ed, depending on the orbital separa- 
tion under consideration. Note that this expression is quadratic in 
e, as expected since the minimum energy corresponds to e = 0. 

Eq. |4T]can be rewritten as the product of two functions. The 
first one caixies the dependence on the orbital separation and goes 
as (X oot'^'i, for a > a^, and oc u)t'-^\ for ot, flc ^ a < fld- The 
second one oc A^ Fn{£,out) / pc {I3 — li/h) depends only on the 
NS structure. A standard calculation provides the value of F„ (^out) 
for different polytropic eos's, while evaluating the integrals in eq. 
1401 is more complex. It first requires to derive the functional form 
of K2{r). The details of this important calculation are deferred to 
our appendix [A] Here we state the main result, 



K.2{r) 



- y{r) - ^VJr) 



3 + y + Ay{r) 



(42) 



where the function y{r) is determined by the solution of the per- 
turbed Poisson's equation for the tidal field (see appendix |AJ. The 
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Table 1. Numerical values of the main structural properties of different 
polytropic eos's, as explained in the text. F„ stands for F„ (§out )■ 



n 


5out 


rj' 


Cl 


K2 


h 


h 


h 





2.449 


0.58 


0.27 


0.75 


198.41 


343.66 


595.23 


1/2 


2.753 


0.43 


0.29 


0.45 


391.99 


817.52 


\1 12.94 


2/3 


2.874 


0.39 


0.30 


0.37 


493.35 


1102.29 


2482.19 


1 


TV 


I/tt 


0.32 


0.26 


785.01 


2023.93 


5299.92 



Table 2. Numerical values of the physical quantities required to calculate 
the free energy in each NS (eqs. l44l and l45t . The three cases refer to different 
choices for the binary components (see text). Different superscripts identify 
the two component stars. Numerical subscripts mean that the corresponding 
quantities are in units of the given power of 10. Units are in c.g.s. 



case/q 


pSs 


Pc,15 


i^d.S 


1^0,3 


p(d) 




A/ 1 


1.27 


1.27 


3.15 


5.31 


~4.2 


~11 


B/0.7 


1.02 


1.27 


2.95 


4.77 


~4.9 


~7.5 


C/1 


1.02 


1.24 


2.91 


5.04 


~2.3 


~2.7 



quantity Aj/(r) — — 3p(r)/p(r) is a measure of the degree of cen- 
tral condensation of the density profile and p{r) = 3A'f (r-)/(47rr'^) 
is the average matter density within r. 

Tab. [T] summarizes our main results for three different poly- 
tropic indices corresponding to plausible NS models. Results for a 
constant-density (n=0) model are also given for reference. 

We are now in the position to quantify the total free energy 
that is pumped in a NS by the tides, during its late inspiral (eg. 1411). 
Since the details of the orbital evolution depend on whether the 
primary is itself a NS or a BH (cfr. sec. |4j, we will discuss the two 
cases separately, beginning with NS-NS systems. 

7.2.1 Free energy in a NS-NS binary 

Let us consider the two stages of their orbital evolution separately. 
During the GW-driven inspiral phase the total energy accumulated 
up to the dynamical radius ad is 



E_ 



(d) 



^°° AEi = (1/2) (^h - ^) j A{h,a^lel^) , 



where we recall that eq^ — {tvluo/uio oj* 
orbital evolution is given by. 



(l + q)i/3- 



We thus obtain, 



E't 



(d) 



34/3 



MX' 



rj4/3 



Stt^F^ 



Fa, and the GW-driven 



(43) 



96 c 

5 



OJd 
LOT 



(44) 



where ^out and F = F„(^out) are given in Tab. [T] Note the pres- 
ence at a very high power of the ratio between the orbital frequency 
at Od and the tidal frequency. This ratio will be in general much 
smaller than unity. Thus, a much larger amount of free energy will 
be accumulated during the subsequent stage of orbital evolution, 
which starts at Od and ends at the contact radius ttc. 



E_ 



(c) 



dE, = (1/2) ( /3 - I 



d (/b,o^o4) (45) 






/3 - -^ 1 ^ ^. 



1-1^ 



In order to derive numerical estimates of Ef in either case 
above, we need to specify the primary mass A/i, the eos and calcu- 
late the dynamical and contact distances. Here, we will concentrate 
on three different types of systems: case A will represent two iden- 
tical NSs with M, = 1.4Af0, i? = 12 Km and n = 1; case B will 
represent a mixed system where the primary is instead a massive 
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Figure 7. As Figure|6]right panel, but for the angular frequency of the wave, 
cjij. In the y — axis we show the fractional distance of ui^, to the orbital 
frequency Wq. 



NS, with the same radius but Af, = 2Mq and n = 1/2; case C 
will include the same massive primary and a massive secondary, 
with n — 2/3. Tab. [2] reports the numerical values of all quanti- 
ties needed in eq. |45] to evaluate the free energy pumped in each 
NS. The first is the available free energy once ad is reached, and 
is expressed in units of 10**^ ergs. The second is the available free 
energy when the contact radius is eventually reached, and it is in 
units of 10**^ ergs. The corresponding energy in the primary can 
be readily obtained, by means of eq. |45] and quantities in the two 
tables. Note however, that its free energy will never exceed that in 
the secondary. 



7.2.2 Free energy in a BH-NS binary 

When the primary is a BH, the ISCO is the larger of the charac- 
teristic radii, if the mass ratio is below some gmax, with 0.3 < 
^max ;$ 0.5, depending on the NS eos (cfr. Fig. 0. This corre- 
sponds to black hole masses larger than a few A/q. This range 
includes the vast majority of BH masses in NS-BH binaries, and 
we therefore focus on this case. 

As for NS-NS systems, the orbital evolution is split into two 
phases, with the GW-driven inspiral being overtaken by a direct 
plunge once aisco is reached. Therefore, the free energy, i5'^'^°, 
available down to aisco will still be given by eq. |44l when we sub- 
stitute ojd with cjisco = i^o(aisc( 



c/aisco\/(l + q)/6. Values 



of this free energy for four different NSs, with a fixed Mi 
black hole companion are reported in tab. |3] 



TMp 



12 Dall'Osso & Rossi 



Table 3. Numerical values of the physical quantities required to calculate 
the free energy pumped in the NS (eqs. l44l and l461 . The three cases refer to 
different choices for the NS polytropic index, while Mbh = 7 Mq is fixed, 
fixing Oist-o — 6.22 X 10® cm (Mbh,©/^)- Numerical subsciipts mean that 
the corresponding quantities are in units of the given power of 10. Units are 
in c.g.s. 



n 


Mns,0 


Pc,15 


Qisco/aT 


'<^T,3 


p(isco) 
^f,45 


f(T) 


1 


1.4 


1.273 


1.42 


3.638 


~ 12.1 


-1.6 


2/3 


1.4 


0.865 


1.36 


3.420 


~3.0 


~1.3 


2/3 


2 


1.236 


1.53 


4.230 


~3.5 


~ 1.4 


1/2 


2 


1.015 


1.48 


4.039 


~ 1.3 


~0.8 



For BH masses < 14Mq, the merger occurs at ax (not Oc), 
which is thus the typical case for stellar mass compact binaries. The 
calculation of the free energy in this plunging phase needs solving 
Sf^ = (1/2) (73 - {Il/h)) /;j^_^ d (/b.oc^o'd) where - con- 
trary to eg. 145 1- the extremes of integration are analytical functions 
of the masses and _R*. We can thus derive an exact solution as a 
function of such parameters, 



E, 



(T) _ Stt^cI F^ M,c- 



Sout L. 

AL_ f O.17 
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l^g^(l + g)VH /3 



h 



1AM, 







0.1 



1- 



where f{uj) = 

The numerical value in the second step is obtained for the 
usual fiducial NS with an n = 1 polytropic eos. Values of the 
free energy for four different NSs and a fixed Mbh ~ 7 Mq are 
reported in the last column of Tab. |3] The total energy at ut is 
~ 10''^ erg. We note that these values depends on the BH mass via 
g*(l + (jf)^. For a 5A/q BH, the available energy would be larger 
by a factor around 4, while for Mi = IOMq it would be reduced 
by a similar factor. 

Following our calculations, (Section 7.2.1 and 7.2.2), we can 
conclude that in both NS-NS and BH-NS systems, the total free 
energy available is ~ 10"'' of the total energy in the tides Eh,T 
(eq.|26t. Therefore, not only the total energy involved in the tidal 
interaction is substantial but also the free energy budget. If tapped 
by additional physical processes, it may have observable conse- 
quences. We devote the following discussion to this possibility. 



8 DISCUSSION AND CONCLUSIONS 

In this paper, we investigate a mechanism of tidal interaction, that 
appears to be most relevant in close compact binary systems. We 
showed that, even in the limit of zero viscosity, a mutual tidal torque 
is maintained successively by the GW-driven orbital evolution and 
by the final dynamical instability before merger. The orbital energy 
transferred to the bulge is quite substantial: 10""^ — 10~^ of the 
secondary binding energy. 

This result is in agreement with Goldreich & Nicholson 
(1989). They consider an inviscid and self-gravitating fluid, subject 
to an external potential, which rotates rigidly at a given frequency. 
They find that, in their case, there is no secular variation of the an- 
gular momentum of fluid particles. Indeed, we showed that, for a 
constant orbital separation, an immediate synchronization occurs 



after which the torque vanishes. The torque is only sustained by the 
orbital shrinking, and it is thus related to secular evolution of the 
total orbital angular momentum 

We then studied the radial structure of the rotating quadrupo- 
lar perturbation, excited within the secondary star by this orbital 
evolution-driven torque. We concluded that, quite generally, it in- 
herits the differential structure of the tidal torque. In the secondary 
frame, the perturbation can be described as a wave excited by the 
motion the the primary. The crest of the wave at different radii ro- 
tates at different rates: while the inner regions promptly follow the 
motion of the primary, the outer regions lag progressively behind. 

This is not a minimum energy configuration, that would cor- 
respond to a uniformly rotating quadrupole. For different binary 
systems, we calculated the free energy associated with this excited 
state, and found that typically ~ 10""^ — 10~' of the total tidal en- 
ergy should be removed from the bulge in order to reach the equi- 
librium state. This free energy reservoir coiTesponds to 10'*^ — 10** 
erg. The most favourable case (10** erg) corresponds to an equal 
mass, double NS system with n = 1. 

A few important considerations lead us to conclude that our 
estimates represent a lower limit on the total free energy. First, we 
adopted the most conservative description of the plunging phase, 
which assumes that the lag angle remains constant during the free 
fall (sec. 17.11 ). In fact, this angle is more likely to increase as the 
orbital separation decreases, thus providing a factor of a few more 
' ^ free energy than we estimated. For instance for our BH-NS systems 

(1 + qfq*f{uj} (46) (see Tab.^, the gain would be a factor (aisco/flT)^ ~ 2.5 - 3.6, if 

the lag angle were to grow linearly with (j;o(a), as suggested by its 
definition (eq. I35l >. Second, our treatment is newtonian and, third, 
linear in the lag angle. Both these features are known to underesti- 
mate the effects of forces and torques, thus the free energy genera- 
tion, close to coalescence, where most of the energy is produced. 

We expect this free energy to be transferred to internal mo- 
tions, via excitation of low-frequency g-modes well before coa- 
lescence and of more effective f-modes, just prior to coalescence, 
when the orbital frequency approaches the fundamental one, a;*. 
We defer to a future work a detailed study of this problem, but it 
is reasonable to expect that the azimuthal motions excited by the 
differentially rotating quadrupolar perturbation will largely inherit 
its differential structure. 

This suggests that a pre-existing magnetic field in the NS core 
might get significantly sheared and amplified. The free energy in 
the excited internal motions grows extremely fast with orbital fre- 
quency and, at the tidal radius, it would correspond to a - mostly 
azimuthal - magnetic field B^ ~ 10^^ G, permeating the whole NS 
volume. As mentioned above, we are probably underestimating the 
free energy content and thus the magnetic field, which nevertheless 
is already as strong as in the most magnetized objects known. In 
fact, this dynamically-driven tidal torque may be able to "recycle" 
an old NS into a magnetar. 

The question is: would the free energy be tapped fast enough 
to build such a strong field component by the time the tidal ra- 
dius is reached? A detailed analysis of this scenario is clearly be- 
yond our scope here, as it would require solving for the fluid mo- 
tions and for their coupling to the magnetic field. However, there 
is an argument that suggests that this may be possible. In a pure 
fi-dynamo the azimuthal field would grow linearly and become as 
strong as the initial poloidal one, before the dynamical radius is 
reached (a > Od). Thus, where most of the free energy is gener- 
ated (a ^ Od), we might expect the linear approximation to break 
down and field growth to proceed in a non linear regime. A growth 
timescale shorter by a factor of 5 — 10 than in the linear regime 



tidal torque in compact binaries 1 3 



Love number radial profile for different polytropes 




^ = It r/R. 

Figure Al. The radial profile of the NS Love number, K2(5), for four dif- 
ferent polytropic eos's, calculated through eg. lA4l 



would be already sufficient to tap the entire free energy reservoir. 
As a matter of fact, when the azimuthal field carries most of the 
magnetic energy, conditions are such that shear instabilities may 
be excited by the differential rotation, (e.g. Balbinski 1985; Luyten 
1990; Shibata, Karino & Eriguchi 2002; Watts, Andersson & Jones 
2005), and the resulting turbulent fluid motions would lead to a 
much faster growth of this field component. 

The above discussion indicates possible measurable conse- 
quences during the final coalescence of a binary system contain- 
ing a NS. If this non-linear phase happens, such a strong magnetic 
field may manifest itself in a magnetar-like flare at the moment the 
NS is tidally disrupted, as it may actually have been observed al- 
ready (cfr. Troja et al. 2010). Alternatively, the magnetic field may 
permeate in a toroidal fashion the NS debris, that will form an ac- 
cretion disc around the remnant of the merger This can favour the 
magnetic production of jets, like it is thought to happen for short 
7-ray bursts. 



APPENDIX A: THE RADIAL PROFILE OF THE TIDAL 
LOVE NUMBER 

The tidal Love number K2 for NSs has been calculated in fully rel- 
ativistic models (Hinderer 2008, D'amour & Nagar 2009, Poisson 
& Binnington 2009, Postnikov et al. 2010), allowing for different 
choices of the equation of state. Building on these works, we sketch 
here a simple derivation of its radial profile within the NS, K2 (r) , in 
the newtonian limit. We start from the relativistic formulation (Hin- 
derer 2008 and references therein), where the spacetime metric is 
written as the sum of a (known) background and a perturbation, in 
axial symmetry. The metric perturbation contains the Laplacian of 
the tidal field of the primary, V'^tpT (cfr. sec.|2]l, and the induced 
quadrupole moment of the secondary, Q22, as leading orders terms. 
Their ratio defines the Love number K2 of the star (Hinderer 2008 
and references therein) 



3 G Q22 

K.2 = -- 



2 Rl V^V'T ■ 



(Al) 



An appropriate choice of the reference frame and gauge conditions, 
reduces the Einstein's equations for the metric perturbation to a 
single, second-order differential equation for the radial function, 
H{r), which describes the tt-component of the metric perturbation. 



Solutions for H{r) are then found separately inside the NS, 
(given a specific equation of state), and in the vacuum outside it. 
The outside metric is matched asympotically to the Schwarzschild 
solution, while a regular behaviour as r — > is imposed to the 
interior solution. The two solutions must eventually be matched 
at the NS surface and this requires that both H{r) and H' (r) be 
continuous across it (cfr. Hinderer 2008, Postnikov et al. 2010). 
This matching determines the coefficient K2, as a function of the 
value of the interior solution at r = i?» . 

By taking the newtonian limit of the relativistic solution, i.e. 
the limit for small stellar compactness, and introducing the vari- 
able y{r) = rH'ir)/ H(r), this procedure eventually recovers the 
classical result (cfr. Brooker & OUe 1955) 



jv J- 

«2 =2 



2-y 
3 + y 



(A2) 



where y here takes the value given by the interior solution for H{r) 
at r = _R, . To find the run of K2 with radius, we need to cut a 
smaller sphere of radius r < R, within the tidally deformed NS, 
and remove the outer ring without changing the matter distribu- 
tion within. Then, we can calculate the Love number of the smaller 
sphere, K2{r), by matching the interior solution with the exterior 
solution at radius r. However, since the matter distribution within 
the smaller sphere is unchanged, it is still determined by the full 
polytropic model. Hence, p{r) 7^ and the matter density has a 
discrete jump at the surface of the sphere, from pij-) to zero. In 
presence of a density discontinuity, the variable y{r) itself will have 
a discrete jump, Ay, which is determined by the new matching con- 
dition at radius r. 

Following Postnikov (2010), it is possible to write this jump 
as 



Ay = 3 



p{r 



p{r~ 



p{r) 



(A3) 



where r and r~ indicate respectively an external point and an 
internal point, infinitesimally close to the surface of the sphere. 
In our case, p{r^) = 0, p{r~) = p{r). Finally, p{r) — 3M{< 
r)/(47rr^) is the average matter density within radius r . Therefore, 
the appropriate function to insert into (ea. lA2L is y(r) + Ay, with 
Ay given by ea. lA3l With this procedure we obtain, 

'2-y{r)-Ay{ry 



i^2{r) 



(A4) 



-> R, and goes to 3/4 (the value 
0. These are the correct limit^ 



3 + y + Ay{r) 

which gives expression ( IA2t for r 
of a constant density sphere) as r - 
for the function K2 {r). 

A generalization of this procedure to the relativistic case is 
conceptually straightforward. 
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will inevitably appear at the boundary of this model. Adding this correction, 
as we have described it here, gives the value of 3/4 also for the calculation 
by Postnikov (2010). 
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